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Abstract. Let V be a separable Hilbert space, possibly infinite dimensional. Let St(p, V) be the Stiefel 
manifold of orthonormal frames of p vectors in V , and let Gr(p, V) be the Grassmann manifold of p dimen- 
sional subspaces of V. We study the distance and the geodesies in these manifolds, by reducing the matter 
i to the finite dimensional case. We then prove that any two points in those manifolds can be connected by 

^ , a minimal geodesic, and characterize the cut locus. 

Resume. Soit V un espace de Hilbert separable, eventuellement de dimension infinie. Soient St(p, V) 
l'ensemble des systemes orthonormes de p vecteurs de V, appele la variete de Stiefel, et Gr(p, V) l'ensemble 
des sous-espaces vectoriels de V de dimension p, appele la variete Grassmannienne. En reduisant le problcmc 
en dimension finie, nous montrons que dans ces espaces il existe des geodesiques minimales entre chaque 
paire de points et nous caracterisons le cut-locus. 



1. Introduction 

Let V be a separable Hilbert space, let p be a positive natural number. We assume that dim(V) > (2p) 
from here on. St(p, V) is the set of orthonormal frames of p vectors in V. Equivalently, we consider 

St(p, V) = {x e L{W, V) : x T o x = ld RP } 

to be the set of all linear isometric immersions of MP into V. Here x T £ L(V,MP) is the transpose with 
respect to the metrics on V and R p , i.e. 

(x T (v),r) Rp = (v,x(r)) v for all v e V, r E R p . 

St(p, V) is a smooth embedded submanifold in V p . The induced Riemannian metric on St(p, V) is (x, y) = 
tr(x T y). St(p, V) is a complete Riemannian manifold with this metric. Gr(p, V) is the manifold of p- 
dimcnsional linear subspaces of V and equals the orbit space St(p, V)/0(p) with respect to 0(p) acting on 
St(p, V) by composition from the right. 

Our interest is due to the fact that St (2, V) with V = L 2 ([Q, 1]) is isometric to the space of planar closed 
curves up to translation and scaling, endowed with a Sobolcv metric of order one. The 0(2)-action on 
St(2, V) corresponds to rotations of the curves. Thus Gr(2, V) with V = L 2 ([Q, 1]) is isometric to the space 
of planar closed curves up to translations, scalings and rotations. See [7], [S], [S] and [B]. Any result that is 
proven about the Stiefel or Grassmannian immediately carries over to the corresponding space of curves. 

2. Critical geodesics 

We will call a curve 7 in a Riemannian manifold a critical geodesic if it is a solution to the equation 
Va t 7 = 0, where V is the covariant derivative. Such 7 is a critical point for the action 3(7, 7) dt. 

Proposition 1 (Critical geodesics in St(p, V) ). Let St(p, V) be endowed with the induced metric from V p . 
Let 7 : [0, 1] — > St(p, V) be a path. Then the geodesic equation is 7 + 7(7 T 7) = 0. Solutions to the geodesic 
equation exist for all time and are given by 

(1) ( 7 We At ,7(<)e A *) = (7(0),7(0))exp/ 



v Id A 

where A = 7(0) T 7(0), S = 7(0) T 7(0), and Id is the p x p identity matrix. 
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For V = W l this has been demonstrated by Edelman et al. [H section 2.2.2J0. Going through their proof 
one sees that the same result holds when V is infinite dimensional. 



Proposition 2. Equation |T]) shows that the subspace ofV spanned by the (2p) columns 0/7(4), 7(f) remains 
in the space spanned by the columns 0/7(0), 7(0) for all t. 

This means that, ifW is the subspace ofV spanned by the columns 0/7(0), 7(0), then we can formulate 
the geodesic equation as an equation in St(2, W). Obviously, dim(W) < 2p. 

This also means that, ifj is a critical geodesic connecting x to y, and the space W spanned by the columns 
of x,y is (2p) dimensional, then, for any t, the columns 0/7(4) and 0/7(4) must be contained in W . 



We denote by d{x, y) the infimum of the length of all paths connecting two points x, y in a Riemannian 
manifold. It does not matter whether the infimum is taken over smooth or absolutely continuous paths0. 
We call a path 7 a minimal geodesic if its length is equal to the distance (£(7(0), 7(1)). Up to a time 
reparamctrization, a minimal geodesic is smooth and is a critical geodesic. We will always silently assume 
that minimal geodesies are parametrized such that they are critical. 

Let (M, g) be a Riemannian manifold, and d be the induced distance. When M is finite dimensional, by 
the celebrated Hopf-Rinow theorem, metric completeness of (M, d) is equivalent to geodesic completeness 
of (M, g), and both imply that any two points x,y £ M can be connected by a minimal geodesic. In infinite 
dimensional manifolds this is not true in general. Indeed, in [T] there is an example of an infinite dimensional 
metrically complete Hilbcrt smooth manifold M and x,y £ M such that there is no critical (and thus no 
minimal) geodesic connecting x to y. A simpler example, due to Grossman [3] (see also sec. VIII. §6 in [3]), 
is an infinite dimensional ellipsoid where the south and north pole can be connected by countably many 
critical geodesies of decreasing length, so that the distance between the poles is not attained by any minimal 
geodesic. 

We will show that, even when V is infinite dimensional, any two points in St(p, V) and Gr(p, V) can be 
connected by a minimal geodesic. 

3.1. Minimal geodesies in the Stiefel manifold. 

Theorem 3. Let V be a Hilbert space. Consider a (2p) dimensional Hilbert space W and an isometric linear 
embedding i : W — » V ■ Then i induces an isometric embedding 



(3) Let x, y £ St(p, W), and a minimal geodesic 7 connecting x to y in St{p, W): then i* 07 is a minimal 
geodesic connecting i*(x) to i*{y) in St(p, V). 

(4) The diameter of St (p, V) is equal to the diameter of St (p, R 2p ). 

(5) Any two points x,y £ St(p, V) can be connected by a minimal geodesic 7. Any minimal geodesic lies 
in St(p,U), where U is a (2p) dimensional subspace of V (dependent on 7J. 

(6) Let x, y £ St(p, V). Then y is in the cut locus of x if and only if there is a (2p) dimensional subspace 
W of V and x,y £ St(p, W) such that x = i*{x), y = i*{y) and i*(y) is in the cut locus of i*(x). 

Note that point (5) in the above theorem implies that minimal geodesies can be numerically computed 
using a finite dimensional algorithm; see Sec. 3.3.4 in [B]. 

We will need two lemmas. 

1. [2] credits a personal communication by R. A. Lippert for the final closed form formula JTJ. 

2. Lemma 6.1 in Chap. VIII in [4] can be used to convert any absolutely continuous path to a shorter and piecewise smooth 



3. Minimal geodesics 




path. 
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Lemma 4. Given x G St(p, V), the set of y £ St(p, V) suc/i that the columns of x, y are linearly independent 
is dense in St(p, V). 

Proof. Let U be the linear space spanned by the columns of x,y; if this space is not (2p) dimensional, then 
let ri, . . . rfc be orthonormal vectors that lie in U ± , with k = 2p — dim([7). Up to reindexing the columns 
of y, we can suppose that the columns x±, . . . , x p , yi, . . . , y p -k are linearly independent. For e > small, we 
then define 

\Vi i = l,...p-k 

yi = \ 

\cos(e)yi + sin(e)rj i = (p - k + 1), . . .p 
It is easy to verify that y G St(p, V) and that the columns of x, y are linearly independent. □ 

Lemma 5. TTie theorem holds when V is a Hilbert space of finite dimension n with n > 2p. 

Proof. We prove point (1). Let us consider the subgroup G = 0{i(W)^ of 0(V) that keeps i(W) fixed. 
Then G acts isometrically on St(p, V) as well, and its fixed point set is z* (St(p, W)) . This proves that 
i*(St(p, W)) is totally geodesic in St(p, V). 

To prove point (2) we first note that since St(p, VF) is isometrically embedded in St(p, V), we have 
d^y^i?/) > dy(i*(x),iH,(2/)). We will show the inverse inequality only for the case when the columns of x 
and y are linearly independent. The general case then follows because the set of y such that the columns of 
x and y are linearly independent is dense in W by lemma U and since distances are Lipschitz continuous. 

Since V is finite dimensional, St(p, V) is compact, so by the Hopf-Rinow Theorem z*(x) and can 
be connected by a minimizing geodesic in St(p, V). The columns of i t {x) and i*(y) together span the (2p) 
dimensional space i(W), so we can apply proposition [5] This allows us to write 7 = i» o 7 for a path 7 in 
St (p,W). Then 

dw^y) < len(7) = len(i« 07)= lcn(7) = d v (i*(x),i*(y)) . 

Point follows from point ^ and the equality 

len(i* 07) = lcn(7) = dvr(^>J/) = dv («*(»), . 

Point follows from point (3). Point ^ follows from the Hopf-Rinow theorem and the discussion in 
Prop. 03 

We now prove point ( 6). By definition, y is in the cut locus of x if and only if there is a geodesic 7 in 
St(p, V) with 7(0) = x, 7(1) = y such that 

supjt : Ien(7| [0tt] ) = d v (7(0), 7(f))} = 1. 

(Recall that we write dy for the distance in St(p,V).) Any such geodesic lies in St(p, W) for some (2p) 
dimensional space W. Letting i : W —> V denote the isometric embedding, we can write 7 = i* o 7 for a 
path 7 in St(p, W). Then one has by point (2) that 

Bup{t:len(7| [0 ,t]) = dw(7(0),7(*))} = l • □ 

We now prove Theorem [3] 

Proof. The proof of points (1), (3), (4), (6) works as in the finite dimensional case. We will now prove point 
(2). We have dw(x,y) > dy U*(x), i*(y)), since St(p, W) is isometrically embedded in St(p, V). It remains 
to show the inverse inequality. 

Consider a smooth path £ connecting i sr {x) to i*(y) in St(p, V). We can find finitely many points = 
to < ti < . . . < ife = 1 such that ^| [ti ,t 4+1 ] is contained inside (the manifold part of) a normal neighborhood. 
So £(£j), £(ii+i) can be connected by a minimal geodesic. By joining all these minimal geodesies we obtain 
a piecewise smooth path 77, with len(r;) < lcn(£). Then by repeated application of proposition [5] there is 
a finite dimensional subspace W of V that contains the columns of r]{t) for t G [0, 1]. When necessary we 
enlarge W such that it also contains i(W). Now the finite dimensional version of this Lemma allows us to 
compare St(p, W) to St(p, W), and we get: 

d w (x,y) = dw-(i*(x),i«(y)) < len(r?) < len(£) . 

Since this holds for arbitrary paths £ connecting ^(x) toi*(y) in St(p, V), we get dw{x,y) < dv(i*(x),i*(y)). 
Point (5) now follows by choosing any linear subspace W containing the columns of x, y. □ 
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3.2. Minimal geodesies in the Grassmann manifold. 

Theorem 6. Thru. [3] remains valid when Stiefels are replaced by Grassmannians. Most importantly, for 
any two points x,y G Gr(p,V), there is a minimal geodesic 7 connecting x to y. The same holds for the 
Grassmannian Gr + (p, V) of oriented p spaces. 

We need a Lemma. 

Lemma 7 (Existence of horizontal paths). For any path x : [0, 1] — > St(p, V) there is a path g : [0, 1] — > 0(p) 
such that the path x(t) o g(t) is horizontal, i.e. normal to the 0(p)-orbits in St(p, V). 

Proof. The tangent space at x to the 0(p)-orbit through x G St(p, V) is 

T x (x.O(p)) = {xz:zE o(p)} = {xz:ze L(W, W), z T + z = 0} . 

Thus y G T x St(p, V) is horizontal if and only if tr(y T xz) = for all antisymmetric z. This is equivalent to 
y T x=0 because y T x is antisymmetric, too. Thus the path xg is horizontal iff 

(d t (xg)) J (xg) = g T x T xg + g T x T xg = 0. 

This can be achieved by letting g be the solution to the ODE g = —x T xg. □ 

Note that the length of x{t) o g(t) is smaller than or equal to the length of x(t), with equality if and only 
if x(t) is already a horizontal path. 
We are now able to prove Thm. [SI 

Proof. St(p, V) is a principal fiber bundle with structure group 0{p) over Gr(p, V) = St(p, V)jO(p). We 
prove the existence of minimizing geodesies connecting any two points in Gr(p, V). Take any point x G 
St(p, V) in the fiber over x. The fiber over y is compact since 0(p) is compact. Therefore d(x, •) attains 
a minimum at some point y in the fiber over y. By theorem [3] there is a minimal geodesic connecting x to 
y. This geodesic is horizontal since otherwise it could be made shorter by making it horizontal. (We use 
lemma [7] here.) By the theory of Ricmannian submersions it projects to a minimal geodesic in Gr(p, V). 

The remaining statements simply follow from Thm. [3] by going to the quotient with respect to the O(p)- 
action. For the case of Gr + (p, V), we use the group SO(p) instead of 0(p). □ 
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